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Abstract. The number of plane partitions contained in a given box was shown by MacMa- 
hon [8] to be given by a simple product formula. By a simple bijection, this formula also 
enumerates lozenge tilings of hexagons of side-lengths a, b, c, a, b, c (in cyclic order) and an¬ 
gles of 120 degrees. We present a generalization in the case 5 = c by giving simple prod¬ 
uct formulas enumerating lozenge tilings of regions obtained from a hexagon of side-lengths 
a,b + k,b, a + k,b,b + k (where k is an arbitrary non-negative integer) and angles of 120 
degrees by removing certain triangular regions along its symmetry axis. 


1. Introduction 

A plane partition is a rectangular array of nonnegative integers with the property that 
all rows and columns are weakly decreasing. A plane partition tt = (7rij)o<i<a,o<i<f) can 
be identihed with its three dimensional diagram = {(i, k) : 0 < k < 7^ij}, and hence 
can be viewed as an order ideal of N^ (an order ideal of a partially ordered set is a subset 
/ such that X & I and y < x imply y E I). 

Replacing each point of by a unit cube centered at it and parallel to the coordinate 
axes, the diagram of tt becomes a stack of unit cubes justihed onto the three coordinate 
planes. By projection on a plane normal to the vector (1,1,1) one obtains that plane 
partitions whose diagrams ht inside an a x 6 x c box are identihed with tilings of a hexagon 
of side-lengths a, 6, c, a, 6, c (in cyclic order) and angles of 120 degrees by unit rhombi with 
angles of 60 and 120 degrees (see [3] or [7] for the details of this bijection); we call such 
rhombi lozenges and such tilings lozenge tilings. 
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Figure 1.1. The hexagon i7(5,5,3); its vertebrae are shaded 
and labeled starting from the indicated horizontal reference line. 

MacMahon showed [8] that the nnmber of plane partitions contained in an a x 6 x c box 
(eqnivalently, the nnmber of lozenge tilings of the corresponding hexagon), with a < b, 
eqnals 

(c + l)(c + 2)^ ■ ■ ■ (c + a)“(c + a + 1)“ ■ ■ ■ (c + 6)“(c + b + 1)““^ ■ ■ ■ (c + 6 + a — 1) 

1 ■ 22 ■ ■ ■ a“ ■ (a + 1)“ ■ ■ ■ 6“ ■ (6 + l)“-i ■ ■ ■ (6 + a - 1) ' 

The elegance of this resnlt strongly invites one to hnd generalizations of it. Two natnral 
directions in which attempts have been made are the following. First, one may consider 
ennmerating order ideals of the prodnct of d hnite chains, with d > 4 (the case d = 4 
was already considered by MacMahon [9]). Nnmerical evidence snggests that there is no 
simple prodnct formnla for the answer to this problem, and very little is known abont it 
(see [11,1 §11]). Second, taking the viewpoint of tilings, one may consider the qnestion 
of connting tilings by nnit rhombi of 2(i-gons with parallel and congrnent opposite sides. 
Again, there seem to be very few resnlts on this in the literatnre for d> 4 (see [4]). 

We generalize MacMahon’s resnlt in the case 6 = c by presenting a general family of 
hexagonal regions with holes on the triangnlar lattice, snch that the nnmber of tilings of 
each member of this family is given by a simple prodnct formnla (examples of regions from 
this family are shown in Fignres 1.2-1.4). 

Consider the tiling of the plane by nnit eqnilateral triangles. Dehne a region to be any 
region of the plane that can be obtained as a hnite nnion of these nnit triangles. 

Consider the hexagonal region i7(a, 6, k) with side-lengths a, b+k^ b, a + k, b, b+k, where 
k is a, nonnegative integer, and snppose H{a, b, k) is drawn snch that its base is the edge 
of length a -\- k (see Fignre 1.1 for an example). Then H{a, 6, k) has a vertical symmetry 
axis £. 
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Figure 1.2. The region 2 , 5 , 7 , 8 , 9 )(5, 6 , 6 ). The vertebrae 

contained in it are shaded, and their labels and the reference 
line for the labeling are shown. 

It is easy to see that H{a, b, k) contains precisely k more nnit triangles pointing npward 
than nnit triangles pointing downward. Since a lozenge always covers one nnit triangle of 
each kind, it follows that i7(a, 6 , k) does not have any tilings for /c 7 ^ 0 . 

This can be hxed by removing some nnit triangles from H{a,b,k), so as to restore 
the balance between np- and down-pointing nnit triangles. One natnral way to do this, 
which was in fact onr motivating example, is to remove an np-pointing triangnlar region 
of side-length k from the center of H{a, b, k) (the instance k = 1 of this is treated in [2] 
and corresponds to a problem posed independently by Propp [10] and Knperberg (private 
commnnication)). At least from the viewpoint of onr proof, it tnrns ont that it is easier 
to ennmerate tilings of the more general regions dehned below. 

We call a triangnlar snbregion of H{a,b,k) which is symmetric with respect to £ a 
window. If the vertex opposite the base of a window is above the base we call the window 
a A-window, otherwise we call it a V-window. A window is called even or odd according 
as its side-length is even or odd. 

The nnit triangles of H{a,b,k) along i can be gronped in pairs forming lozenges — 
we call them rhombic vertebrae — except possibly for those snch triangles tonching the 
bonndary of H{a,b,k), which we call triangular vertebrae (instances of both kinds of 
vertebrae occnr in Fignre I.I; the vertebrae are indicated by a shading). The dehnitions 
below involve the concept of labeling the vertebrae starting from a reference line — a line 
of the triangnlar grid perpendicnlar to i: by this we mean that we label the vertebrae that 
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Figure 1.3. The region i 7 ( 2 , 4 ), ( 3 , 5 ) (7, 8 , 3). The vertebrae 
contained in it, their labels and the reference line are shown. 


fnlly he on each side of the reference line starting with 1 for the closest vertebra, 2 for the 
next, and so on; see Fignre 1.1 (in case the reference line is the base of H{a, 6 , k) and the 
vertebra tonching it is triangnlar, we skip this vertebra and start by labeling the next one 
by 1 ). 

Snppose k is even. Remove from H (a, 6 , k) some nnmber of even A-windows of total size 
k (the total size of a collection of windows is the snm of the side-lengths of the windows). 
Since we removed only even windows, each vertebra is either fnlly removed or fnlly left in 
place. Label the vertebrae of i7(a, b, k) starting from its base, and snppose the vertebrae 
in the leftover region have, in increasing order, labels (see Fignre 1.2 for an 

example). Let 1 = (/i,..., In)- The region obtained from H{a, 6 , k) by removing windows 
as above may have forced lozenges (e.g., if some of the windows tonch). Remove all snch 
forced lozenges. Denote the leftover region by i7i(a, b, k). 

If k is odd, start, from top to bottom, by removing a nnmber of even A-windows from 
H{a, 6 , /c), then remove a single odd window (which can be either a A- or a V-window) and 
continne downwards by removing some even V-windows (the set of removed even windows 
above the odd window, or the set of removed windows below it, or possibly both, may 
be empty). In this process, obey the restriction that the total size of the A-windows is k 
larger than the total size of the V-windows. Choose the base of the single odd window as 
onr reference line, and label the vertebrae on both sides of it as indicated above. Snppose 
the vertebrae below the reference line in the leftover region have, in increasing order, labels 
/i,..., Im: and the ones above it have, in increasing order, labels gi,..., (note that one 
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Figure 1.4. The region i^(i,3,4),(i,4)(®’ 1)’ 
with the vertebrae it contains. 


or both of these lists may be empty). Let 1 = (/i,..., q = (gi,..., Remove any 
forced lozenges, and denote the leftover region by iLi,q(a, 6, k) or Lfi^q(a, b, k), according as 
the removed odd window was a A- or a V-window (Fignres 1.3 and 1.4 show an example 
of each kind). 


The regions iLi(a, 6, k), Hi qia, b, k), and b, k) are the regions whose nnmber of 

tilings we will determine, thns generalizing the case 6 = c of MacMahon’s resnlt (which is 
obtained for /c = 0). 


Before giving the formnlas expressing the nnmber of tilings of these regions, we need to 
dehne two families of polynomials, both indexed by pairs of lists of positive integers. 

For m, n > 0 dehne monic polynomials Bm,nix) and Bm,n{x) by the formnlas 

5’ 


= 2—™(—i)/2(a: + n + l)^{x + n + 2)^ 
{x + 2){x + 3)^ ■ ■ ■ {x + n — l)‘^{x + n) 



2n + 1\ 
2 ) 


n 


{x + i). 


Y\{2x + n + i + 2) 


BmA^) = 2-™-(-+i)/2(a: + m + 1)^ 

{x + l)(a; + 2)^ ■ ■ ■ {x + m — l)‘^{x + m) 


n+i—1 


X + 


X + 


X + 


2m — 3 


X + 


2 m — 1 


n 


{x + i). 


Li {x + i + l/2)n fj 


W{2x + m + i + l)m+i, 


( 1 . 1 ) 


( 1 . 2 ) 


where {a)k '■= a(a + 1)... (a + k — 1) is the shifted factorial (in the middle two lines of 
(1.1) and (1.2), the bases grow by one from each factor to the next, while the exponents 
increase by one midway throngh and then decrease by one to the end). 

Next, given (possibly empty) lists of strictly increasing positive integers 1 = (/i,..., Im) 
and q = (gi,..., q^), dehne constants ci q and ci^q by setting 


/ n — m\ 

Cl,q = 21 2 

m ^ n ^ 

Y\.l<i<j<Ah ^i)Y\.l<i<j<nA A 

11 (21.)! 11 (2„ - 1) 1 


Cl,q = 2l 2 )- 

mTJ ^ TT 

'Wl<i<j<Ah ~ ni<i<j<n(^J ~ A 

IJ (21. - 1)! n (2,.), 



(1.3) 

(1.4) 

Let A = A(l) and fi = /u(q) be the partitions having parts /m, —ni, Im-i —m + 1,..., /i — 1 
and Qn — n, Qn-i — n + 1 ,..., gi — 1 , respectively. 

We are now ready to dehne onr polynomials PiAA PiAA- They are given by the 
following formnlas: 


B\,q{x Ai) — Cl qS.fj2,n(^) 

]^(a: — h{c) + m + l)(a; + h{c) + n + 1) ]^(a: — h{c) + n + 2)(a: + h{c) + m) 
ceA ce/^ (1.5) 

B\,q{x Ai) — Cl qS.fj2,n(^) 

]^(a: — h{c) + m + l)(a: + h{c) + n) ]^(a: — h{c) + n + l)(a 


c M- m 


c€A 


cEl^i 
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( 1 . 6 ) 



where for a cell c = (i, j) of a partition diagram, h{c) is defined to he i + j (in all formnlas 
above empty prodncts are considered eqnal to 1 ). 

Denote by 1*-*^ the list obtained from 1 by omitting its i-th element, and by 1 — 1 the list 
obtained from 1 by snbtracting 1 from each of its elements (except in case /i = 1 , when 
1 — 1 is dehned to consist of /2 — 1 , ■ ■ ■ , /m, — !)• 

Onr generalization of MacMahon’s ennmeration of plane partitions contained in a sym¬ 
metric box is given in Theorem 1.1 below. Its proof is dednced in Section 3 from Propo¬ 
sition 2 . 1 , which interprets the above-dehned polynomials Pq,i(a:) and .Pq,i(a;) as tiling 
ennmerators of two families of lattice regions dehned in Section 2. Proposition 2.1 and 
Theorem 1.1 are the main resnlts of this paper. 

Let M{R) denote the nnmber of lozenge tilings of the region R (the “M” stands for 
“matchings,” as lozenge tilings can be identihed with perfect matchings on the dnal graph). 

Theorem 1.1. The number of lozenge tilings of the regions Hi{a,b,k), k), and 

ffi qidj b, k) is given by the following formulas: 

(a) For k even, we have 


2-^M{H^{a,h,k)) 


-P 0 ,i((® -hk — 2 )/ 2 ).Pi_i^ 0 (a/ 2 ), 
-P 0 ,i((o^ + k — 2 )/ 2 )Pi_ij(a/ 2 ), 

P 0 i(m) ((a k — l)/2).Pi 0 ((a — l)/2). 


a even, li = 1 , 
a even, li > 1, (1-7) 

a odd. 


( 6 ) For k odd, we have 


2-—M(77i,q(a,6,fc)) 


Pl,q((a + /c - l)/2)Pq^ 

l(^) (a/ 2 ), a even, 
A,q(») ((« + ^)/ 2 )Tq,i((a - l)/ 2 ), a odd, 


( 1 . 8 ) 


and 


2-m-nM(H,,^(a,b,k)) 


Pi,q((a -|-/c — l)/2)Pq i(,7i) (a/ 2 ), a even, 
-Pi,q(") ((® + ^)/2).Pq,i((a - l)/2), a odd. 


(1.9) 


Remark 1.2. For k ^ 0, the bijection between tilings of onr regions and stacks of nnit 
cnbes breaks down. Indeed, trying to lift snch a tiling to three dimensions one is lead to 
an “impossible” constrnction, similar to Penrose’s impossible triangle (see Fignre 1.5). 

In Section 3 we show how to rednce the statement of Theorem 1.1 to the problem 
of hnding the tiling generating f un ction of certain simply-connected regions, dehned in 
Section 2. In Sections 4 and 5 we prove, nsing an indnctive argnment, that these tiling 
generating fnnctions are given by the above-dehned polynomials P and P. The idea of the 
proof is very simple: nsing combinatorial considerations, we dednce recnrrence relations 
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Figure 1.5. A tiling of i7(i)^(i^4^5)(5,5,3). 


for the tiling generating fnnctions of onr regions, and we show that the polynomials P and 
P satisfy the same recnrrences. This approach, however, raises the qnestion of how conld 
one gness that the tiling generating fnnctions of onr regions are given by the intricate 
formnlas (1.1)~(1.6). We explain in Section 6 how this can be done in a conceptnal way, 
nsing some compnter assistance. 

2. Two families of regions 

We now dehne two families of regions that will tnrn ont to be closely related to the 
polynomials P and P given by (1.5) and (1.6). 

Consider the triangnlar lattice drawn so that one of the directions of the lattice lines 
is horizontal. From each point of the lattice, there are six possible steps one can take on 
the lattice. We call the fonr non-horizontal steps, according to their approximate cardinal 
direction, northeast, northwest, sonthwest and sontheast. 

Let O be a point of the triangnlar lattice and consider the inhnite lattice paths Pu 
and Pd that start from O and take alternate steps northeast and northwest, respectively 
sonthwest and sontheast (these paths are shown in dotted lines in Fignre 2.1(a)). We 
say that two consecntive steps of Pu, the hrst going northeast and the second northwest, 
form a bump] similarly, a step sontheast followed by a step sonthwest on Pd form a bnmp. 
Label the bnmps of Pu and Pd, starting with the ones closest to O, consecntively by 
1 , 2 ,.... Select arbitrary hnite snbsets of the bnmps of Pd and Pu, and snppose their 
labels, in increasing order, form (possibly empty) lists 1 = {li,... ,1^) and q = [qi,... ,qu), 
respectively. Assnme that at least one of 1 and q is non-empty. 

Starting from these selected bnmps, we dehne the “right bonndary” of a region i?i q(a:) 
as follows (here a: is a non-negative integer on which a lower bonnd will be imposed below). 
Call the nppermost and lowermost points on a bnmp the top and bottom of that bnmp, 
respectively. For each selected bnmp on Pu, draw a horizontal ray to the left of its top, 
and a ray going sonthwest from its bottom. For each selected bnmp on Pd, draw a ray 
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Figure 2.1(a). i?(2,4,5),(2,4)(2)- 


going left from its bottom and a ray going northwest from its top. 

By following these rays, we join the selected bnmps on in a connected piece: follow 
the horizontal ray from the top of bnmp qi to its intersection point with the ray starting 
from the bottom b of bnmp Qi+i, and then nse the latter ray to reach b. The selected 
bnmps on Pd are connected into a single piece in a similar way. 

To join the above two pieces together, we draw a horizontal ray left of O. We follow 
the ray going sonthwest from the bottom of bnmp qi of Pu nntil its intersection with the 
horizontal ray originating at O. Then we follow this latter ray nntil it hits the ray going 
northwest from the top a of the /i-th bnmp of Pd, which we then track nntil we reach 
a. We have now inclnded all selected bnmps in a connected path Q joining the top A of 
the last selected bnmp on P^ to the bottom B of the last selected bnmp on Pd (in case 
1 is empty, B is taken to be the intersection of the ray going sonthwest from the bottom 
of the hrst selected bnmp on P^ with the horizontal throngh O; for q = 0 , A is chosed 
analogonsly; see Fignres 2.1(b) and 2.1(c)). It is easy to see that this path has exactly 
2lm — m + n + 1 steps going sontheast (the one exception is the case 1 = 0 , when it has 
2lm — m + n = n steps). 

To close the bonndary of i?i,q(a:), we continne from B by taking x steps left (in case 
1 is empty, we move left nntil we are x -\- 1 nnits from O). Then we head northwest for 
2lm — m + n + 1 steps (one fewer when 1 = 0 ), after which we tnrn and go northeast nntil 
we hit the horizontal ray going left from A (see Fignre 2.1(a)). Finally, take x-\-{lm — Tn) — 
{qn — n) + 1 steps east on this ray, nntil we reach A. It is not hard to see that the length 
of the sonthwest portion of this bonndary mnst be chosen as indicated (i.e., to match the 
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number of southeast steps on Q) in order for the enclosed region to admit lozenge tilings 
(this follows for example by encoding tilings as families of non-intersecting lattice paths, 
as described at the beginning of Section 4). 

In general we will allow the lozenges in a tiling to be weighted. The weight of a tiling 
is the product of the weights on its lozenges. The sum of weights of all tilings of a region 
R is denoted M{R) and is called the tiling generating function of R. 

We dehne to be the region enclosed by the above-dehned boundary, requiring in 

addition that the tile positions htting in the selected bumps on Pu are weighted by 1/2 (i.e., 
if one of these positions is occupied by a lozenge in some tiling, that lozenge gets weight 
1/2 — this will be indicated in our hgures by a shaded oval placed on such tile positions; 
all other tiles get weight 1). Since the top and bottom sides have non-negative length, x 
must satisfy x > max{ 0 , — Im — — ^} (and x > max{— 1 , — n — 1 } = qn—n — 1^ 

for 1 = 0). Dehne i ?0 0 ( 2 :) = 0 for all x. 

Our second family of regions, denoted .Ri q(a:), is dehned almost identically (see Figures 
2.2(a)-(c)). The only difference is that, in the construction of the path Q that connects the 
selected bumps, instead of the horizontal ray starting from O, we use the one starting from 
the lattice point (5, one step southwest of O (and, in case 1 is empty, when dehning the 
bottom side of .R 0 q( 2 :), we move left until we are x units away from O). Asa consequence, 
this path has now one fewer southeast steps than before, and hence we dehne the southwest 
boundary of .Ri,q( 2 :) to have length 2lm — m + n. Also, the top side of the enclosed 
region has now length x + (Im — m) — {qn — n). The parameter x satishes therefore 
X > max{0, qn ~ Im — P (this time the condition on x remains the same for 1 = 0). 
Again, we dehne .R 0 0 ( 2 :) = 0 for all x. 

The close connection between the above-dehned regions and the polynomials P and P 
given by (1.5) and (1.6) is expressed by the following result. 

Proposition 2.1. For all lists 1 and q of strictly increasing positive integers we have 
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Figure 2.2(a). i?(2,4,5),(2,4)(3)- 



Figure 2.2(b). i?0,(2,4)(5)- 


Figure 2.2(c). 4,5), 0(3). 


M(i?i,q(a:)) = Pi,q(a:) 

M(.Ri,q(a:)) = A,q(ai), 


for all X for which the regions on the left hand side of these equalities are defined. 




Figure 3.1. 


In the next section we show how the above Proposition implies Theorem 1.1. The proof 
of Proposition 2.1 is presented in Sections 4 and 5. 

3. Reduction to simply-connected regions 

One nsefnl way to approach certain tiling ennmeration problems (lozenge tilings in 
particnlar) is to identify them with families of non-intersecting lattice paths, and then use 
the Gessel-Viennot theorem (see e.g. [12, Theorem 1.2] or [5]) to express the number we 
seek as a determinant. While this could be carried out directly for the regions on the left 
hand sides of (1.7)~(1.9) that we are concerned with, it turns out that it will be more 
effective to reduce hrst our problem to the case of certain simply connected subregions, 
for the number of tilings of which we can easily obtain recurrences. This reduction can be 
achieved using the Factorization Theorem for perfect matchings presented in [1, Theorem 
1 . 2 ], 

Let R be a region of the triangular lattice, and suppose R has a vertical symmetry axis £ 
(an example is shown in Figure 3.1). Let Pi and Pr be the two zig-zag lattice paths taking 
alternate steps northeast and northwest and touching i from left and right, respectively. 
Group the unit triangles of R crossed by £ in sequences of contiguous triangles. These 
are separated by sequences of triangles of the lattice not contained in R-, we call these 
sequences gaps. 

We cut R into two subregions and R~ as follows. Gonsider the paths Pi and Pr. 
Start at top by following Pr to the hrst gap. After each gap, keep following the path along 
which we reached it or switch to the other path, according as the gap contains an even or 
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an odd number of unit triangles. Continue this until we reach the bottom of R. Finally, 
assign weight 1/2 to all tile positions in R that have two sides contained in our cutting 
path. Dehne R'^ and R~ to be the regions obtained to the left and right of our cutting 
path, respectively, and having the above-mentioned lozenge positions weighted by 1/2. 

Since we are interested in counting lozenge tilings, we may assume that R has an even 
number of unit triangles crossed by £ (indeed, if this number was odd, the symmetry of R 
would imply that R contains an odd number of unit triangles, so it would have no lozenge 
tilings); dehne the width of i?, denoted w(i?), to be half this number. 

Then the Factorization Theorem for perfect matchings of [1, Theorem 1.2] implies 

M{R) =2^^^'^M{R+)M{R-). (3.1) 

Indeed, lozenge tilings of R can be identihed with perfect matchings of the dual graph 
G of R, i.e., the graph whose vertices are the unit triangles of R and whose edges connect 
precisely those pairs of triangles that share an edge. Since R is symmetric with respect 
to £, so is G, and all conditions in the hypothesis of [1, Theorem 1.2] are easily seen to be 
met. We obtain that M(G) = 2^ M(G'’") M(G“), where G^ and G~ are certain precisely 
dehned subgraphs, and M(G) denotes the matching generating function of G. However, it 
is easily checked that G^ and G~ are exactly the dual graphs of the regions R'^ and R ~, 
respectively, thus proving (2.1). 

Proof of Theorem 1.1. All equalities in the statement of Theorem 1.1 follow directly by 
applying (2.1) to the regions on the left hand side of (1.7)-(1.9) and using Proposition 2.1. 
More precisely, in all cases, the absolute value of the exponent in the power of 2 in (1.7)- 
(1.9) turns out to be the width of the corresponding region, and, by Proposition 2.1, 
the two factors in the product on the right hand side turn out to be equal to the tiling 
generating functions of the two regions that arise by applying the Factorization Theorem 
to our regions. 

(a) Take R to be the region i7i(a, b, k), and apply (2.1). Since all the gaps are even, the 
cutting path determining 7?+ and R~ lies fully on the right side of i (this is illustrated in 
Figures 3.2(a)-(c)). Suppose a is even. Since a -|- /c is also even, this cutting path starts 
at the midpoint of the top side of R and ends at the midpoint of its base (see Figures 
3.2(a) and (b)). Consequently, all unit triangles of R crossed by i can be paired to form 
rhombic vertebrae. By our dehnition, these are precisely the /i-th,..., /m-th vertebrae of 
the hexagon H enclosed by the outer boundary of R. Therefore, the region 77+ to the left 
of the cutting path is precisely 7^0 i( (a + k — 2)/2). 

As for the region R~ to the right of the cutting path, suppose hrst that /i > 1 (see 
Figure 3.2(a)). Then it is easy to check that the region obtained from R~ after removing 
the forced lozenges is precisely .. ;^_i)(a/2), rotated by 180 degrees. On the other 

hand, if /i = 1, because of the different pattern of forced lozenges in R~ (extra forced 
lozenges on the bottom, as shown in Figure 3.2(b)), the subregion left after removing 
the forced lozenges is isomorphic to 77(;2-i,...,Zrn-i)(®/2)- Since the width of R is m, this 
proves, by Proposition 2.1, the hrst two equalities in (1.7). 

For a odd we proceed similarly. Now the last labeled vertebra of 77 is a triangular 
vertebra, and it remains present for all choices of removing even A-windows from 77 
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Figure 3.2(a). 
-f^( 2 , 3 , 4 , 6 , 7 )( 6 , 5, 4) reduces to 
-R0,(2,3,4,6,7)(4) and i?(l, 2 , 3 , 5 , 6 ), 0 ( 3 ). 



Figure 3.2(b). Figure 3.2(c). 

Hu, 2 , 4 , 6 , 7 )( 6 , 5, 4) reduces to i 7 (i^ 3 ^ 4 ^ 6 j)( 5 , 5, 4) reduces to 

-R 0 ,( 1 , 2 , 4 , 6 , 7 )( 4 ) and .R(1,3,5,6),0(3)- -R 0 ,( 1 , 3 , 4 , 6 )( 4 ) and .R( 1 , 3 , 4 , 6 , 7 ), 0 ( 2 )- 


(see Figure 3.2(c)). By definition, its label is Im- This fact and an examination of the 
pattern of forced lozenges in R'^ and R~ explains why, at the indices of the two factors 
on the right hand side of the third equality in (1.7), the list 1 appears once with its last 
entry omitted and once in full. Indeed, because of forced lozenges at the top of i?"*", the 
triangular vertebra labeled Im is not relevant for the region obtained from 7?+ by removing 
the forced lozenges. This leftover region is easily seen to be exactly 7^0 ((a + /c — l)/2). 

Furthermore, the region obtained from R~ by removing the forced lozenges is isomorphic 
to .Rij((a — l)/ 2 ), irrespective of the value of h (the difference in this regard from the 
case a even is explained by the fact that the pattern of forced lozenges in R~ depends on 
/i for a even, but does not for a odd). Since the width of R is m, by Proposition 2.1 this 
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Figure 3.3(a). 
jf^( 2 , 4 ),( 3 , 4 )( 6 , 8 , 1 ) reduces to 
-R(2,4),(3,4)(8) and i?(3,4),(2)(3)- 


Figure 3.3(b). 
;f^( 2 , 4 ),( 3 , 4 )('^’ 8 , 1 ) reduces to 
-R(2,4),(3)(4) and i?(3,4),(2,4)(8)- 



^( 2 , 3 ),(i, 3 , 5 )( 6 , 8 , 1 ) reduces to 
-R(2,3),(1,3,5)(3) and i?(i,3,5),(2)(3). 



^(2),(i, 3,5,6)(7,8, 1) reduces to 
-R(2),(1,3,5)(4) and i?(i,3,5,6),(2)(3). 


completes the proof of (1.7). 

(b) To prove (1.8), take R to be the region i7i^q(a, b, k) and apply (2.1). The resulting 
cutting path stays on the right of I until it arrives to the odd A-window, then it switches 
to the left of £ and stays there until it reaches the bottom of R (this is illustrated in 
Figures 3.3(a) and (b)). 

Recall that our reference line for labeling the vertebrae of H is now the base of the odd 
removed window; denote it by 77. Suppose a is even. Then the last vertebra below 77 is in 
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this case a triangular vertebra, which remains present for all choices of removed windows 
(see Figure 3.3(a)). The rest of the unit triangles of R crossed by I can be paired to form 
rhombic vertebrae. By dehnition, the region R contains precisely the vertebrae labeled 
li,... ,lm below R, and those labeled qi,... ,qn above R. It is easy to check that the region 
obtained from R'^ after removing the forced lozenges is precisely .Ri^q((a +/c — l)/2), while 
the one obtained from R~ after the same procedure is congruent to i?q pm) (a/2). Since 
the width of i? is m + n, we obtain by Proposition 2.1 the hrst equality of (1.8). 

If a is odd, a similar analysis shows that R'^ is congruent to .Rj q(n) ((a +/c)/2), and R~ 
to Rqy((a — l)/2) (see Figure 3.3(b) for an illustration). Since the width of R is again 
m + n. Proposition 2.1 implies the second equality of (1.8). 

The regions .Rpq(a,6,/c) are treated in a perfectly similar way (this is illustrated in 
Figures 3.4(a) and (b)). Again, the regions R"*" and R~ obtained by applying (2.1) turn 
out to be members of the families of regions dehned in Section 2, and applying Proposition 
2.1 one arrives at the formulas (1.9). □ 

Remark 3.1. Part (a) of Theorem 1.1 (i.e., the case k even) also follows from Lemma 2.2 
of [6]. Indeed, as shown in the next section, it is easy to express the tiling generating 
functions of the regions R"*" and R~ arising by applying the Factorization Theorem to 
Ri(a, 6,/c) as determinants of Gessel-Viennot matrices (see (4.1)). It is also not hard to 
check that the entries of these matrices are of the form required by [6, Lemma 2.2], which 
therefore provides explicit product formulas for the value of these determinants. 

For k odd on the other hand, the corresponding Gessel-Viennot matrices have a more 
complicated structure and their determinants cannot be computed by [6, Lemma 2.2]. 

4. Recurrences for M(Rpq(a:)) and M(Rpq (x)) 

Let 1 = (/i,..., Im) and q = (gi,..., q^) be two (possibly empty) lists of strictly in¬ 
creasing positive integers, and let T be a tiling of Rpq(a:) (Figure 4.1 shows an example 
for 1 = (2,3), q = (2,4) and a: = 3; a white center on a tile indicates that it is weighted 
by 1/2). This tiling divides the southwestern boundary of Ri,q(a:) into 2lm — m + n + 1 
unit segments. Starting from these segments, by following the tiles of T, one obtains 
paths of rhombi that end on the southwest-facing unit segments on the right boundary 
of Ri,q(a:) (these paths are shaded dark in Figure 4.1). It is clear that these paths are 
non-intersecting, and it is not hard to see that they completely determine T. Since along 
the paths of rhombi the only possible steps are northeast and east, it is clear that they 
can be identihed with lattice paths on the square lattice taking steps north and east. 

We will need the following special case of the Gessel-Viennot theorem for non-intersec¬ 
ting lattice paths (see [5]). All our lattice paths will be on the grid graph Z^, directed so 
that its edges point east and north. We allow the edges of 7? to be weighted, and dehne 
the weight of a lattice path to be the product of the weights on its steps. The weight of 
an V-tuple of lattice paths is the product of the individual weights of its members. The 
generating function of a set of V-tuples of lattice paths is the sum of the weights of its 
elements. 

An V-tuple of starting points u = (wi, ... ,un) is said to be compatible with an V-tuple 
of ending points v = (ui,..., n^v) if for alH < j and k < I, every lattice path from Ui to 
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Figure 4.1. Encoding a tiling as paths of rhombi. 
vi intersects each lattice path from Uj to Vk- 

Theorem 4.1 (Gessel-Viennot). Let u and v be compatible N -tuples of starting and 
ending points on the graph 7? oriented as above. Then the generating function for N-tuples 
of non-intersecting lattice paths joining the starting points to the ending points is 

det ((ciij)<iv) 5 

where aij is the generating function for lattice paths from Ui to Vj. 

We will hnd it convenient to view the paths of rhombi directly as lattice paths on Z^, 
thns bypassing the “extra steps” of bijecting them with lattice paths on a lattice of rhombi 
with angles of 60 and 120 degrees, and then deforming this to the sqnare lattice. In this 
context, the “points” of onr lattice C are the sonthwest-facing edges of the triangnlar 
lattice — we call them segments —, and the “lines” of C are seqnences of adjacent rhombi 
extending horizontally or in the sonthwest-northeast direction — we call these rows and 
columns, respectively. The rhombi that make np these rows and colnmns are called edges 
of £. 

To coordinatize onr lattice, we choose the x-axis to be the bottommost row, and the 
y-axis to be the leftmost colnmn, intersecting i?i^q(a:). Let N = 2/^, —m + n + 1, and label 
the starting and ending segments of the lattice paths encoding onr tiling T, from bottom 
to top, by wi,..., Wat and ni,..., n^v, respectively. 

Weight by 1/2 the edges of C corresponding to tile positions weighted by 1/2 in i?i q(a:). 
Weight all other edges of C inside i?i^q(a:) by 1, and the edges ontside i?i^q(a:) by 0. By 
Theorem 4.1 and the above-mentioned bijection between tilings and lattice paths, we 
obtain that 

M(i?i,q(a:)) = det ((aij)i<ij<iv), 

where aij is the generating fnnction for lattice paths from Ui to Vj. 
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(4.1) 



Suppose m <n. Then expanding this determinant along the last row turns out to give 
us a recurrence relation for M(i?i q(a:)). 

Lemma 4.2. For m < n, we have 

n 

M(i?i,q(a:)) = M(i?i,q(.) {x)) + (-1)^5™ M(^i,q(a:)), (4.2) 

k=i 

where Smn is the Kronecker symbol and 

^ / X + X Im F Qk \ / . „N 

^ \2q'fc + m —ny 2\2qk + m — n — lJ 

Proof. Let A be the matrix in (4.1). Expanding the determinant in (4.1) along the last 
row we obtain 


N-l 

M(i?i,q(x)) = {-l)’^aN,N-k det{A^’^-^), (4.4) 

k=0 

where A'^i is the matrix obtained from A by deleting row i and column j. 

We claim that at most the hrst n + 1 terms in the sum (4.4) are nonzero. Indeed, recall 
that ON^N-k is the generating function for lattice paths on C going from the segment un 
to the segment VN-k- It follows from our construction of i?i,q(a:) that there is a unique 
segment s on its right boundary contained in the row immediately below O, and s is 
facing southeast. The ^-coordinates of the segments un and s work out to be, in our 
coordinatization, 2/^ — m -|- n and 2/^, respectively. Since our paths on C take steps 
northeast and east, and m < n, we see that a path starting ai un will end, on the right 
boundary of i?i,q(a;), at some segment above s, or, possibly, in case m = n, at s itself (see 
Figures 4.2(a), 4.3(a) and 4.4(a)). Since there are precisely n segments above s on the 
right boundary of i?i^q(a:) (one corresponding to each selected bump on P^), this proves 
our claim (and also implies s — vn-u)- 

Furthermore, it is easy to determine the exact value of aN,N-k for k < n. It follows from 
the previous paragraph that aN,N-n = 0 unless m — n, in which case aN,N-n = 1 (this 
accounts for the Kronecker symbol in (4.2)). For 1 < k < n, ON^N-n+k is the generating 
function of lattice paths going from un to VN-n+k- In our system of coordinates it turns 
out that we have 


Un = (0, 2lm — m + n) 

VN-n+k = (x + Ijn - Qk - m + n + l,2qk + 2lm), for 1 < k < n. 

Partitioning the lattice paths from un to VN-n+k in two classes according to the direc¬ 
tion of their last step we obtain that aN,N-n+k = Ck, with Ck given by (4.3). Replacing 
the upper summation limit in (4.4) by n and then replacing the summation index k by 
n — k we obtain from (4.4) that 
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Figure 4.2(a). 


Figure 4.2(b). 




M(J?i,q(a:)) = det(A'^’'^--+") + (-l)"5^n det(A'^’'^-"). (4.5) 

k=i 

To complete the proof of the Lemma, we show that the determinants on the right hand 
side of (4.5) are eqnal to the corresponding tiling generating fnnctions on the right hand 
side of (4.2). 

More precisely, we claim that, for 1 < /c < n, we have det(A^’^“"^+^) = M(i?i ^(fe) (x)). 
Indeed, by Theorem 4.1 it follows that det(A^’^“"^"'“^) is eqnal to the generating f un ction 
of {N — l)-tnples of non-intersecting lattice paths starting at ui,... ,un-i and ending 
at vi,... ,VN-n+k-i,VN-n+k+i, ■ ■ ■ ,vn- However, by the bijection between tilings and 
lattice paths, this is precisely the generating f un ction for tilings of q(fe)(a:) (for 1 = (2,3), 
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Figure 4.4(a). 


Figure 4.4(b). 


q = (2,4), this is illustrated in Figures 4.2 and 4.3; the starting and ending segments of 
the paths on the left are omitted, leading to the regions on the right). 

Similarly, one sees that = M(.Ri^q(a:)). Indeed, is the gener¬ 

ating function for {N — l)-tuples of non-intersecting lattice paths going from ui,..., un-i 
to ui,..., VN-n-i, VN-n+i, ■ • •, Vjsf, respectively. In this case the corresponding region has 
the same sets of selected bumps as i?i^q(a;), but its boundary around O is changed so that 
it becomes precisely .Ri^q(a:)) (see Figures 4.4(a) and (b)). This completes the proof of the 
Lemma. □ 

The above result gives a recurrence relation for M(i?i q(a:)) when m < n. Note that this 
assumption was indeed necessary: for m > n there are additional non-zero terms in the 
sum (4.4), and the regions they can be associated to are not part of our families R and R. 
To deal with the case m > n we use the “symmetry” of i?i,q(a:): we encode its tilings by 
non-intersecting lattice paths starting on its northwestern side and taking steps southeast 
and east. This leads us to the following result. 

Lemma 4.3. For m > n, we have 


m — l 

M(i?i,q(a:)) = + + (4.6) 

k=i 


where 


Dk 


l'x + lm + lk-'ITl + n + l\ 
\ 2lk-m + n + l ) 


(4,7) 


Proof. We proceed in the same fashion as in the proof of Lemma 4.2. First, we change 
the dehnition of the lattice C so that its points (called again segments) are the unit edges 
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Figure 4.5(a). 


Figure 4.5(b). 


of the triangular lattice facing northwest, and the lines (still called rows and columns) are 
sequences of adjacent unit rhombi extending horizontally and in the northwest-southeast 
direction. We choose the x- and y-axis to be the topmost row and the leftmost column 
intersecting Ri^qix), respectively. Since the northwestern side of i?i.q(a:) has length 2qn -t- 
m — n, a tiling of i?Lq(a;) is encoded by this many non-intersecting lattice paths in C. Let 
N = 2qn + m — n and label the starting and ending segments of these lattice paths, from 
top to bottom, by wi,..., un and ui,..., respectively. 

Formula (4.1) holds also with this new dehnitions of £, Ui and Vi. Expanding the 
determinant on its right hand side along the last row we obtain 


N-l 

M(i?i,q(x)) = {-lfaN,N-k (4.8) 

fc =0 

where the entries of A are now dehned in terms of the new UiS and ViS. 

We claim that all terms in this sum except for the hrst m are 0. Indeed, denote by s the 
unique segment in the right boundary of Ri^q{x) contained in the row immediately above 
O. The ^-coordinates of un and s in T are 2qn + m — n — 1 and 2qn — 1, respectively. 
Since m > n, and our paths take steps southeast and east, it follows that the path starting 
at 'Utv ends at a segment below s on the right boundary of i?i^q(a:) (this is illustrated in 
Figures 4.5(a) and 4.6(a) for 1 = (2, 3, 5), q = (2, 4) and x — 2). There are exactly m such 
segments, one for each selected bump on P^. So aATfc = 0 unless k G {N — m -t- 1,..., N}, 
which proves our claim. 

For 1 < /c < m, the exact value of aN^N-m+k follows from the coordinates of 'Utv and 
VN-m+k- Since we have 
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Figure 4.6(a). 


Figure 4.6(b). 


un = (0, 2qn + m — n — 1) 

VN-rn+k = {X + lm- h, 2/fc + 2^^), foT 1 < k < Ul, 

we obtain that ttN^N-m+k = where is given by (4.7). Therefore, replacing the 
npper snmmation limit in (4.8) by m — 1 and then replacing the snmmation index k by 
m — k, (4.8) becomes 


m— 1 

M(i?i,q(rr)) = det(yl^’^-"^+^) + Dm det(yl^^). (4.9) 

k=i 

To complete the proof of the Lemma, we show that the determinants on the right hand 
side of (4.9) are eqnal to the corresponding tiling generating fnnctions on the right hand 
side of (4.6). 

We claim that, for 1 < /c < m — 1, we have det(A^’^“™+^) = M(i?i(fe) ^(a; — 1)). 
Indeed, by Theorem 4.1, det(A^’^“”^"'“^) is eqnal to the tiling generating fnnction for 
(A^ — l)-tnples of non-intersecting lattice paths in C starting at wi,... ,un-i and ending at 
ni,..., VN-m+k-i, VN-m+k+i, ■ • ■, vn- By the bijection between tilings and lattice paths, 
this is precisely the generating fnnction for tilings of i^pfc) q(a: — 1). (The change in the 
parameter a; is a new ingredient that was not present in the case m < n; this is dne to the 
change of onr coordinate system in the present case. Fignre 4.5 illnstrates this for k = 2.) 

Similarly, one sees that det{A^’^) = ,^{x + Im — Im-i ~1))- Indeed, det{A^’^) 

is the generating fnnction for (N — l)-tnples of non-intersecting lattice paths going from 
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ui, ... ,un-i to vi,... ,vn-i, respectively. The corresponding region is the portion of 
^l,q( x) consisting of the nnit triangles in and above the row containing vn-i, and to the 
right of the colnmn containing un (see Fignre 4.6). This is readily seen to be exactly 
i?i(m) q(a: + Im — Im-i “ 1)), thns completing the proof. □ 

We will also need the analogs of the above two Lemmas for the regions .Ri^q(a:). 
Lemma 4.4. Form < n, we have 


M(^,,q(a:)) = 5^(-l)'^-"CfcM(4q(.)(a:)), 
k=i 

where 


f X + Im -\- Qk + + 

\2qk + m — n + 1/ 2 \2qk + m — nj 


(4.10) 


(4.11) 


Proof. Let C be as in the proof of Lemma 4.2. The sonthwestern side of .Ri,q((r) has 
length N = 2lm — m + n. Encode the tilings of .Ri,q((r) as Wtnples of non-intersecting 
lattice paths on C. Order their starting and ending points from bottom to top and express 
their nnmber as a determinant by Theorem 4.1. Using the same reasoning as in the proof 
of Lemma 4.2, it follows that only the last n entries in the last row of this determinant are 
nonzero. These are easily seen to be eqnal to the nnmbers Ck given by (4.11). Expand 
this determinant along the last row. Using argnments similar to those in the proof of 
Lemma 4.2, the order N — 1 minors in the determinant expansion can be interpreted as 
tiling generating fnnctions of regions from the .R-family, and one obtains (4.10). □ 

Lemma 4.5. For m >n, we have 


ra — 1 

M(R,,q((r)) = {-l)^-^Dk M(Ri(.)_q((r - 1)) + M(Ri(^),q((r + /^ - - 1)) 

k=i 

+ (-l)"^h^^ M(Ri,q((r - 1)), (4.12) 

where 6mn is the Kronecker symbol and 


Dk 


X -\-1 

m + lk - m + n 
2lk — m + n 


(4.13) 


Proof. Choose C to be as in the proof of Lemma 4.3. The northwestern side of .Ri.q(a:) 
has length N = 2qn + m — n + 1. Encode the tilings of Ri,q(a:) as Wtnples of non¬ 
intersecting lattice paths on C. Order their starting and ending points from bottom to 
top and express their nnmber as a determinant by Theorem 4.1. By argnments similar to 
those in the proof of Lemma 4.3, it follows that only the last m -|- dmn terms in the last 
row of this determinant are nonzero. More precisely, it tnrns ont that the (m -|- l)-st to 
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last entry in the last row is dmm and the snbseqnent entries are given by (4.13). Expand 
the determinant along the last row. In analogy to the sitnation in the proof of Lemma 4.3, 
it tnrns ont that the minors of order — 1 in this expansion are eqnal to tiling generating 
fnnctions of regions belonging to the R- and /^-families. Precise identihcation of these 
regions leads to (4.12). (Note that, except for the last term on the right hand side of 
(4.12), the shift in the parameter x is precisely the same as in the case of Lemma 4.3. The 
last term does not have a correspondent in (4.6) becanse there we didn’t allow m = n.) □ 

Finally, to snccessfnlly carry ont the indnctive evalnation of onr tiling generating fnnc¬ 
tions, we need to pay attention to their behavionr when 1, q and x are on the “bonndary” 
(i.e., 1 or q is empty, or x takes on the minimnm valne for which onr regions are dehned). 
By dehnition, we take M(i ?0 0 ( 0 :)) = 1 and M(.R 0 0 ( 2 :)) == 1, for all x (empty regions have 
precisely one tiling). 

In the recnrrences (4.2), (4.6), (4.10) and (4.12), the valne of the argnment in the terms 
on the right hand side is at least as large as one less the valne of the argnment on the 
left hand side. Therefore, if at least one of 1 and q is nonempty, these recnrrences can be 
applied as long as the argnment x on their left hand side is not minimnm possible — i.e., 
for (4.2) and (4.6), if x is strictly larger than max{0, Qn — Im — n + m — 1} (max{—1, — 

n — 1} = qn — n — in case 1 = 0), and, for (4.10) and (4.12), if x is strictly larger than 
max{ 0 , qn — Im — + w}- 

When X does take on its minimnm possible valne, we can express the tiling generating 
fnnctions of onr regions as shown in the following two Lemmas. 

Lemma 4.6. (a) Suppose Im — ^ ^ Qn — n. Then for 1 7 ^ 0 we have 

M(i?i,q(0)) = M(i?iM,q(/^ - - 1)). (4.14) 

(6) Let Im — 'ITT' R ^ Qn — 'n. Then we have 

M(i?i,q(Q„ - Im - n + m - 1)) = M(i?i q(„) {q^-l^-n + m-l)). (4.15) 

Proof. Recall that, for 1 7 ^ 0, the lengths of the base and top side of Ri^q(a:) are x and 
X -|- (fm — m) — {qn — n) -t- 1, respectively. The assnmption in (a) implies therefore that 
the base is at most as long as the top side. Thns, by setting a: = 0 we obtain a valid 
region, which has all the lozenges along its sonthwestern side forced (see Fignre 4.7(a) for 
an illnstration). The leftover region is readily identihed to be Ri(m) q{lm — Im-i — !)• 

Consider now part (b). Onr assnmption implies that the top side is shorter than or 
eqnal to the base. Indeed, for 1 7 ^ 0 this follows from the above paragraph. For 1 = 0, the 
top side has length x — {q^ — n) + 1, and the base has length x — qi + 2. Therefore, it 
follows that the top side is at most as long as the base also in this case. 

By setting x = qn — Im — 'n + 'm — 1 the top side shrinks to a point and we obtain 
a valid region whose lozenges along the northwestern side are forced (see Fignre 4.7(b)). 
The leftover region is easily seen to be Ri^q(n) {qn — Im — tx + m — T). Since the top forced 
tile has weight 1/2 and the others have weight 1, we obtain (4.15). □ 




Figure 4.7(a). Figure 4.7(b). 


Lemma 4.7. (a) Suppose Im — ^ Qn — n. Then for 1 7 ^ 0 we have 

M(i?l,q(0)) = M(.RiM,q(/^ - - 1)). (4.16) 

( 6 ) Suppose Im — TTb < Qn — n. Then for q 7 ^ 0 we have 

M^Ri^diqn - Im - n + m)) = ^ M(.Ri q(„) {q^ - Im - n + m)). (4.17) 

Proof. For 1 7 ^ 0, the lengths of the base and top side of .Ri^q(a:) are x and x + {Im — 
'^) ~ (in ~ respectively. Therefore, the assnmption in (a) implies that the base is at 
most as long as the top side. Thns, by setting x = 0 we obtain a valid region, which has 
all the lozenges along its sonthwestern side forced (see Fignre 4.8(a)). It is not hard to see 
that the leftover region is precisely — Im-i — !)■ 

For part (b), note that onr assnmption implies that the top side is at most as long as 
the base. Indeed, we have already seen this for 1 7 ^ 0. For 1 = 0, the top side and the base 
have lengths x — {qn — n) and x — qi + 1, respectively. Therefore, the top side is as least 
as long as the base in this case as well. 

By setting x = q^ — Im — n + m, the top side shrinks to a point and we obtain a valid 
region whose lozenges along the northwestern side are forced (see Fignre 4.8(b)). The 
leftover region is readily seen to be .Ri q(Ti) {qn — Im — n + m). Since the top forced tile has 
weight 1/2 and all others have weight 1, this proves (4.17). □ 

5. Proof of Proposition 2.1 

In this section we show that the polynomials Pi,q{x) and Pi^q(a;) satisfy the same re- 
cnrrences as M(Ri.q(a:)) and M(.Ri.q(a:)). Then we dednce Proposition 2.1 by indnction 
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Figure 4.8(a). 


Figure 4.8(b). 


on Im + n + X. 

Rewrite formulas (1.5) and (1.6) as follows: 


m li — 1 

P\,ciix) = ClqBm,nix + Im - m) Yl + Im - j) (x + Im - m + U + j + 2) 

i=l j=i 
n Qi — l 

n n (3^+- TO+n - j+ 1 ) (x++j+ 1 ) 

j=i j=i ( 5 . 1 ) 

m li — 1 

Pl,Ciix) = CmBm,n{x + /m “ to) (x + - j) (x + - m + U + j + 1) 

i=l j=i 
n Qi — l 

n n (31 +- TO + n - j)(a: +/^n + J + 1 ). 

i=i j=i (5.2) 

Lemma 5.1. Form < n, we have 

n 

Pl.^{x) = 5^(-l)"-"C'fcPi,q(.) (x) + {-iThmuPxA^), (5.3) 

k=l 

where 5mn is the Kronecker symbol and Ck is given by (4.3). 

Proof. The outline of our proof is the following. Dividing both sides of (5.3) by 

CnPi,q{x), the terms on the right hand side turn out to become ratios of linear poly¬ 

nomials in x, all having the same numerator. Dividing through further by this common 
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numerator, the identity to be checked becomes the partial fraction decomposition of the 
rational function on the left hand side. Since the denominator of the latter is a product 
of distinct linear polynomials in x, this identity can be readily verihed. 

The expression on the right hand side of (4.3) giving Ck can be rewritten as 

^ {2x + 2lm-m + n + 2){x + lm-qk-m + n + 2)2q^+m-n-i 

2{2qk + m-n)l ' ^ 

The left hand side. We have 


Pi,ci _ Crq _ 2 {2qn + m - n) ! _ 

CnP\ ci(n) Cj q(?i) (2x T 2/^ fn T n T 2)(x T Im qn T n T 2^2qn+m—n—l 

n™ 1 +im - j)ix+im - m+n+j + 2 ) 

]\T=1 rijC /{x + Im - + Im - m + n + j + 1) 

nr=i njLV(^ +i^+j + i){x+i^-m+n-j+i) 

nr=i^ {x + Im + j + i){x + Im - m + n - j) 

Pm,nix T l^n TiT) 

Pm,n—lix T Im 


/ n — m\ , 1 

2( 2 )—m+l 


i2qn + m — n)l 


2 
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n l TT 1 Y[l<i<j<mi^j ni<j<j<n(^J T) 


1 


n — 1 


1 


0(2'.) lip*-!)! 


n 


l<i<j<m 


(ii - h) n 


l<i<j<n-l Vyj 


-liQj - Qi 


nr=inS 


n—1 


-p-r X + lm—Tn + n + li + 1 -p-r X + lm—IXl + U — i + 1 

^\x + lm-m + n + i + l ^\x + lm-m + n-qi + l 
1=1 1 = 1 

Qn 1 

Y[ (x + Im + j + 'i-)ix + Im - m + n - j + 1) 


J=n 

1 

^m-\-n 


m — 1 


(2x “h ‘^Ijn ^ Tl 


x + lm + n — i + l 


x + lm + ri-i + 1 / 2 ' 

1=0 


where we used the (easily checked) fact that 


Pm,nix') _ 

Pm,n-lix) ~ 2^+^ 


m — 1 


^i2x + m + n + 2)m+n ^ 


x + m + n — i + 1 


a: + m + n — t + 1/2 

i=0 


(5.5) 


(6.6) 
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After cancelling ont common factors in the nnmerator and denominator of (5.5), and 
then rearranging the remaining factors, we obtain 


P\,ci _ i‘2qn + m - n) ! {qn - Qi) 

— 1) ! 11^=1 (^i + Qn) 

(2x + 2lm -m + n + 2)n-ni I\T=i{x + Im + h - m + n 1) 
{2x + 2lm - m + n + 2) nr=/ [x + Im - Qi - m + n + 1) 

The terms in the sum on the right hand side. Using (5.1) and (5.2) we obtain 


P'kPl.qW 

P'nPl^qin) 


Ci^q(fc) {2qn + m 


Ci,q(-) {^Qk + m 


Kzl n 


9i+l- 

3=i 


n 


9i —1 
j=i 


- n)! (x + - Qfc - m + n + 2)2q^+m-n-i 

- n) \ {x + Ira - Qn - m + n + 2)2g„+m-n-l 

^{X + Im+ j + l){x + lra-m + n- j) 

{X + Irn + j + 1){X + Im - m + n - j) 


{2qn + m- n)l 
{2qk + m- n)\ 

m ^ r 

iit^ n 


ni<i<j<Tn Y\.l<i<j<rL, i,j^k [Qj - Qi) 




nr=in;=i,,^.(/.+?.) 


n —1 




Y\.l<i<j<mih ^i)Y\.l<i<j<n-l^Qj Qi) 


UT=iUUih + Q,) 


^ (2b)! 5^(2,.-!)! 

{x -|- Irn Qk 'kkl n 2)2q^.-\-m—n—l 
{x A Im Qn Xn n 2)2q^-\-m—n—l 

n—1 gi+i — 1 

n n {x + Irn + j + ^)ix + Irn — 1X1 + n — j). 

i=k j=qi 

By cancelling ont common factors and regronping one obtains 


(5.8) 


CkPi,q(k) _ {2qk - 1 )! {2qn + m-n)l k + qk 

CnPi,qM (2<?n “ 1 ) ! {^Qk + m - n) ! b + qn 

]Ti=li^kiQn - Qi) _ x + lm-Qn-m + n + 1 

itiZl(Qk - Qi) U7=k+1 iqi-qk)^ + ^rn-Qk-m + n + l 

The last term on the right hand side. Since for m 7 ^ n this term is 0, assnme m = n. 
We have 
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-1 


^l,q ^ Cl,q _ 2 {2qn) ’ _ 

^n-Pl,q(") q(n) {2x + 2/^ + 2) (x + /^ — Qn + 2)2g„ 

nr=i rijC/ +/n - j) {x+in+j+1 ) 
nr=l 11^=/ ( 2 ^ + - j) (^ + + J + 1) 
nr=i - j){x +/„+j+1) 

nr=i^ + L-j){x+in+j + i) 

Bn,nix In Tl') 

Bn,n—lix In 


2 ( 2 )-".+! 

2( 2 )“"■ 


n 


(29,.)! 




n 


1 Yll<i<j<nih ^i)Yll<i<:j<ni^j 


nr=in;=i 


n—l 


iiT^n 


ni<i<j<n-l(9j 9i) 


(2!.)! (29i - 1)! 

fJri.”! I 


U7=lUUih + q,) 


nj=n (^ + — j)(a: + /n + J + 1) 

{2x + 2ln + 2,)ix + In — Qn + 2 ) 2 q„-l 
ix In fl 1)2715 


(5,10) 


where we have used the equality 


Bn,nix^ 
Bn,n—l (^) 


[x + 1)2715 


(5.11) 


which can be readily checked. Simplifying and rearranging in (5.10) leads to 


1^1,q ll ' ' ' In D[j=l i^n Qi) X -\- In Qn “1“ 1 

^7i-Pl,q(") Ql'''(ln-1 nr=l(^i “I" X + In + 1 


(5.12) 


By (5.7), (5.9) and (5.12), it follows that the statement of the Lemma is equivalent to 
the following rational function identity: 
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{2qn + m - n) ! nr=/ (gn - Qi) 

{2qr,-l)l Utiih + Qn) 

{2x + 2lm -m + n + 2)n-m + Im + h - m + n + 1) 

{2x + 2lm -m + n + 2) Y\a=i {x + Im - Qi - m + n + 1) 


n 


= B-i) 


X + Ir 


(Xk _ 

-Qfc-m+n+l 


+ 


X + In - qn + ^ 


(5.13) 


where at and f3 are given by 


{2qk - 1) ! {2qn + m - n) ! YlT=ii^i + Qfc) 117=1 - Qi) 

{2qn - 1) ! {2qk +m-n)l UT=iih + Qn) - qt) ]T^Zl+l{q^ - Qk) 

h---ln nr=/(9n-9i) 

qi---qn-i Yl7=i(^i + (In) ' 


However, (5.13) can be regarded as the partial fraction decomposition of the rational 
fnnction on its left hand side. Since the denominator is factored in distinct linear factors, 
the valnes of the constants ak and 13 that make (5.13) trne can be immediately obtained 
by snitable specializations of x. It is easy to see that these valnes are precisely the ones 
in the above two eqnalities. □ 

Lemma 5.2. Let Dk, Ck and Dk be given by (4.7), (4.11) and (4.13). 

(a) For m > n, we have 


771—1 

Pi,ci{x) = ^ (-l)""“^IlA;Pi(fc)^q(a: - 1) + + /^ - Im-i - !)• (5.14) 

k=i 

(6) For m < n, we have 


hAx) = (5.15) 

k=l 

(c) For m> n, we have 


m—1 

-Pl,q(3:^) = ^ DkPl(k) q{x — 1) + DmPl(m) q{x + Im — lm-1 ~ 1) 

k=l 

+ i-irSmnPl,ci{x-l). 
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Proof. The identities (5.14)-(5.16) can be proved nsing the same approach as in the 
proof of Lemma 5.1. For the hrst two, we divide throngh by the last term on the right 
hand side, and for the third by the second to last. In each of the three cases, the terms 
on the right hand sides become ratios of linear polynomials, all with the same nnmerator. 
Dividing throngh by this common nnmerator, the identity to be proved becomes the partial 
fraction decomposition of the rational fnnction on the left hand side. In all three cases, 
the denominator of this rational fnnction is a prodnct of distinct linear factors, and thns, 
as in the proof of Lemma 5.1, the identity is easily checked by snitable specializations of 

X. 

This can be carried ont in a straightforward way. The shifts in the argnments in 
(5.14)-(5.16) are so that the division mentioned above resnlts in great simplihcation of the 
involved expressions. The only difference that needs to be pointed ont is that we need the 
following relations, analogons to (5.6) and (5.11), involving the monic polynomials Bm,n 
and 


^ 1 + 'm + n){x + m + n + l){2x + m + n + 2)m+n-i 


Bm,n(.X^ _ 1 

S^,n-l(^) 


n-1 . 

n x + m + I 
X + m + i + l/2 

i=0 

■{x + m + n){2x + m + n + l)m+n 


x + m + n — i — 1 


Bm,n{^^ 
Bnri—l,n{x^ 2 


m —1 

n X -\- m -\- n — i — 1/2 

i=0 

^ n — 1 

^(2a: + m + n + l)m+n JJ 


i=0 


X + m + i + 1 
X + m + i + 1/2 


Bn,n{x - 1) _ {x)2n+l 


Bn—l,n(^x'j X fl 

These are all easily verihed directly nsing the dehning formnlas for Bm,n and Bm,n- D 
Lemma 5.3. We have 


Pi,q{qn -lm-n + m-1)^ {qn - Im - n + m - 1) 

Phqiqn -Im-n + m) = ^ A,q(M (in “ ^ + m). 


-Pl,q(0) — Tj(7ti) q(/j7j Iyti — I 1) 

-Pl,q(0) ~ Tj(7ti) q(/j 7 ^ Iyti — I !)■ 
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and for 1 ^ 0 we have 



Proof. For m < n, from (5.7) and (5.4) one obtains a simple expression for the ratio 
It is readily seen that for x — qn — lm. — 'n+m — l this expression becomes 
1/2, thns proving the hrst eqnality in the statement of the Lemma in this case. For m > n, 
one obtains from (5.5) an expression for Pi^q(a:)/(C'nPi^q(n)(x)) similar to (5.7). Again, it 
follows easily that the ratio Pi,q(a:)/Pi q(r.)(a;) specializes to 1/2 for a: = qn — lm—n + m — l. 

The remaining three eqnalities are proved similarly, jnst as the proof of Lemma 5.2 is 
similar to that of Lemma 5.1. □ 

Proof of Proposition 2.1. We proceed by indnction on Im + n + x. For i?i^q(a:), the 
minimnm possible value of /^ + n + a: is — 1. In this case 1 and q are necessarily empty, 
and (2.1) is true by our dehnitions. For .Ri^q(a:), the minimum value of /^ + n + a: is 0. 
Again, this implies 1 = q = 0, and (2.2) follows from our dehnitions. 

Let N > 0 and suppose that (2.1) and (2.2) are true for all 1, q and x with Im+n+x < N. 
Consider 1, q and x so that Im + n + x = N. We claim that (2.2) holds also for these 
choices of 1, q and x. 

Indeed, assume hrst that x is not equal to its minimum possible value. Suppose m > n. 
By (4.12), M(.Ri.q(a:)) can be expressed in terms of tiling generating functions of regions 
from the R- and (in case m = n) i?-family. By the induction hypothesis, (2.1) and (2.2) 
hold for all the regions on the right hand side of (4.12). Then Lemma 5.2(c) implies that 
(2.2) holds in this case. For m < n, (2.2) follows similarly, from (4.10), the induction 
hypothesis and Lemma 5.2(b). 

Assume now that x is equal to its minimum possible value. Then x equals the value 
of the argument on the left hand side of either (4.16) or (4.17), and we can express 
M(.Ri q)(a;) in terms of the tiling generating function of a region from the .R-family for 
which the induction hypothesis applies. Therefore, Lemma 5.3 implies that (2.2) holds for 
^l,q( x) as well. 

We now show that (2.1) holds for Im + n + x = N. Assume that x is not equal to 
its minimum possible value. Suppose m < n, and apply Lemma 4.2. By the induction 
hypothesis, (2.1) holds for the regions appearing in the terms of the sum on the right hand 
side of (4.2). Furthermore, we have seen above that (2.2) holds for .Ri^q(a:). Therefore, by 
Lemma 5.1 we obtain that (2.1) is true for Ri^qix), if m < n. The case m > n follows 
analogously, using Lemma 4.3, the induction hypothesis and Lemma 5.2(a). 

Finally, if x takes on its minimum possible value, (2.1) follows by (4.14) and (4.15), 
the induction hypothesis and Lemma 5.3. This completes the proof of Proposition 2.1 by 
induction. □ 


6. The guessing of M(Ri q(a:)) and M(Ri q(a:)) 

In this section we describe how we arrived at conjecturing that the tiling generating 
functions for the regions RLq(a:) and .RLq(a:)) are given by the polynomials P\^q{x) and 
A,q(a:). 

We discuss here the case of Ri,q(a;). The region Ri q(a:) can be treated similarly. 

It is easy to write down explicitly the entries of the matrix in (4.1). This gives an easy 
way to compute the polynomials M(.Ri q(a:)) by computer, for specihc choices of 1 and q. 
At once, one is striked by the fact that all these polynomials seem to factor in the form 
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M(i?i,q(a:)) = Ci,qFi,q(a:), (6.1) 

where ci^q is a constant and Fi q(a:) is the prodnct of linear factors of the form (x + t), 
where t is either an integer or a half-integer. 

Moreover, it is easy to conjectnre the behavionr of Fi q(a:) nnder incrementing a single 
element of 1 or q. Specihcally, after some experimenting one soon comes np with the gness 
that 


Fiifc)^q(a;) 

,q(^ ~ 

i^i,q(a:) 

Fj q|fe) (x) 

^l,q(^) 


{x - Ik + lm){x + Ik + Im - m + n + 1), 

- == x{x + 2lm — m + n + 1) 

{x + Qk + Im + ^){x - Qk + - m + n), 


for 1 < /c < m 


for 1 < /c < n, 


( 6 . 2 ) 

(6.3) 

(6.4) 


where is the list obtained from 1 by increasing its k-th element by 1 (so for k < m, 
ll^^ is dehned only if Ik+i — Ik >2). These formnlas provide ns with an explicit gness for 
Fi^q, in terms of F[^] [^](a:), where [m] is the list consisting of 1,..., m. Fnrthermore, some 
experimentation easily leads one to conjectnre that F[m],[n]{x) = Bm,n{x) (where Bm,n{x) 
is given by (1.2)). 

It is considerably more difhcnlt to gness a formnla for the constants ci^q jnst based on 
working ont examples (this is not snrprising, as factorization of integers gives signihcantly 
less information than factorization of polynomials). However, it tnrns ont that we do not 
need to gness the valne of these constants: snbstitnting the conjectnred expression for Fi q 
in (6.1) and nsing then Lemma 4.7, we obtain (conjectnred) recnrrences satished by the 

ci,q’s. 

More precisely, snppose Im — rn > Qn — n. We may then nse (4.16) to relate ci^q to 
ci( m) q. Carrying this ont one arrives at 


Cl,q 

Cp-m.) q 


n"7('™-*.) 


■jm—n —1 


(2/m — 1) ! Y\i=l{^-rn + <ii) 


On the other hand, if Im — xn < Qn — n, we dednce similarly from (4.17) that 


(6.5) 


0,q _ pn-m-l ^ 11^=1 i^ri Qi) 

Cl,q(») {‘^Qn) '■ UZMn + k) ' 


( 6 . 6 ) 


We can then repeat this procednre for 1 and q replaced by l^'”^ and q or 1 and q^"^\ 
depending on which of (4.16) and (4.17) applies. At each iteration we have to look at the 
initial segments of the lists 1 and q that are left over, and compare the difference between 
their largest entry and their nnmber of entries, before deciding which of (6.5) or (6.6) 
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to apply. However, it is not hard to see that no matter in which order we apply these 
recnrrences, they lead to an expression of the form 


Cl,q = a'" ’’ n 


n 


1 n 


l<i<j<Tn 


(h-um 




{<lj - Qi) 


{2k - 1)! {2qk I nr=i n;=i(/.+ 

where e(l, q) is some integer depending only on 1 and q. Then after some experimenting 
one is readily lead to gness that e(l, q) = —m, and we obtain (1.4). It is remarkable 

how the “weak” conjectnres (6.1) and (6.2)-(6.4) lead ns, by the “forcing” argnment in 
Lemma 4.7, to the precise (conjectnred) formnla for ci,q. 

A similar analysis can be done for M(i?i^q(a:)), and one arrives at formnlas (1.5), (1.1) 
and (1.3). 


Acknowledgments. David Wilson’s program Vaxmacs for connting perfect matchings 
was very nsefnl in hnding the families of regions i?i^q(a:) and .Ri^q(a:). 
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